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ABSTRACT 

This report develops a simple model to describe the low altitude shear associated 
with a microburst (or small downburat) in. the atmosphere. Such microbursts are thought 
to be the cause of convective radial flows near the ground which represent potentially 
hazardous conditions to aircraft during take-off and during approach and landing. Cloned 
form solutions are presented for the equations of mass, momentum and thermal energy to 
give the spreading rate of the microburst and the velocity and temperature decay due to 
turbulence. An analysis of the impact of the microburst and its radial spreading over the 
ground is given and illustrative examples are provided. The solution for the downward and 
outward wind velocities are given in tenn3 of simple dimensionless algebraic parameters 
which can be computed readily in realtime during a piloted simulation of an aircraft flying 
through a microburst. 
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I. Introduction 

This paper presents a simple model to describe the low altitude shear associated with 
a microburat (or small downburst) of lower temperature air ia the atmosphere. Such 
microburst3 are thought to be the cause of coavectiye radial flows near the ground which 
represent potentially hazardous conditions to aircraft during takeoff and during approach 
and landing. 

A typical microburst has been described 1 as having high downward velocities, perhaps 
up to a maximum of G0-7rt ft/sec, over a region of 1,C00 to 10,000 feet diameter, and lasting 
many minutes in duration. The microburst impinges on the ground and spreads radially, 
causing local horizontal wind velocities comparable to those in the downward flow, and i3 
accompanied by intense shear profiles near the ground. An aircraft flying through such a 
phenomena would experience initially a headwind, and subsequently a tailwind as it passed 
through the microburat. 

During the psol several years a number of attempts have been made to describe the 
velocity and temperature fields associated with downburst phenomena. The most simple 
of these 2 ia a correlation of observed velocity profiles which provides a scmiempirical, easily 
calculated, algebraic expression that is suitable for use in flight simulation. A physically 
more realistic model 3 , also suitable for use in flight simulation, takes account of some of the 
fluid dynamics of the flow by representing the microburst by a distribution of singularities 
(sources, vortices and doublets) but doe3 not take into account the temperature Geld 
and its associated buoyancy effects, nor the effect of turbulent diffusion. More elaborate 
models 4 ’ 5 *® take into account those effects but require substantial computational time to 
provide numerical results and are generally not suitable for real-time simulation. 

The work presented here ia an attempt to develop a model which is suitable for use with 
real time flight simulators, permitting adjustment of parametric constraints to fit observed 
measurements while retaining the essential physical aspects of buoyancy and turbulent 
diffusion. 
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2. Analysis 

The approach used in the present analysis assumes that the microburst is a steady axi- 
ally symmetric downward flow issuing from a point (ie, a virtual origin) in the atmosphere 
at distance h 0 above the ground. It is acted upon by turbulent diffusion and by buoyancy 
associated with a temperature difference (Too~T 0 ) between the center of the downflow and 
the ambient atmosphere. In the close vicinity of the ground the flow is turned in a radial 
direction and forms a ground jet. No account is taken of the effects of ground roughness 
in the present analysis (ie, a smooth ground is assumed). Figures 1-3 show the general 
schematic for the flow and the geometry of the downflow and subsequent radial outflow 
near the ground. 

2.1 The Downflow 

The differential equations for conservation of mass, momentum and thermal energy in 
the presence of gravity are written for an incompressible fluid, a3 
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where and c e are eddy diffusivitie3 associated with the velocity and thermal fields. 
In the foregoing equations temperature ~nd density changes are assumed to be small 
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The integral forma of equations (2) may be found by integration radially across the 
flow, ie, 
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= 0 


( 5 ) 


since the light hand side vanishes at both limits. 

The second term in equation (5) may be written, using equation (1), as 




so that equation (5) becomes 


d r°° 
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( 0 ) 


Equation (6) i3 the integral form of the downward momentum equaiion and describes 
the rate of increase of downward momentum that results from the (nega ive) buoyancy of 
the colder air within the downflow. 

Similarly, the integral form of equation (3) may be found as 
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or, using equation (1), 
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equation (7) becomes: 



constant 


( 8 ) 


Equation (8) is the integral form of the equation of thermal convection in the absence 
of thermal sources. 

It in assumed that the dovvnburst originates in a localised region of the atmosphere 
which can be approximated ao a virtual point source so that all quantities depend on a 
similarity variable. More specifically solutions are considered of the form 


“ =t “(i)" /(f) ' * “ ^ ( 0 “ h({) w 

where f = y/c, and Ui and 0 t arc reference values of u and 6 at the point s — Si measured 
downward along the axis from a virtual origin s = o located at an altitude h 0 . 
Substitution of these expressions into equations (G) and (8) gives 
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Equations (10) and (1 1) yield for the exponents m and n 
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Substitution of these values into equation (10) gives a relation between u, 0, and Zj ie 


u| _3 Jf h(f)fdf , . 

Og*i 4 J?PkMs K * 

It is well known that for turbulent jet flows the profiles /(f) and h(f) are Gaussian in 
character, ie they can be approximated by 


m = 

h(s) 


where a and b are constants. 

Substitution of these profiles into equation (13) gives a relation between a anc b, ie 


( 15 ) 

(The constant b i3 known from experimental results to be of the order of 100) Equation 
(15) indicates that the velocity profile is influenced by the thermal conditions thiough the 
parameter E. 

The constants a and b also established the relationship between the half width of the 
microburst (where the velocity is equal to half the value at the center of the jet.) 
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Thus, setting 
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>-(¥)' 

and the diameter of the microburst (within which the velocity exceeds one h.. f of the value 
at the center) is 

dm -* ,-(iSaa)*. 

or 
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Similarly, the thermal diameter of the nJcroburat is 







To summarise, the solution for the downflow in the microburst can be wiitten in terms 
of reference values Ui and occurring at a location jj 
The velocity and temperature distributions are given by 


and 


where 
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These equations are readily calculated in real-time to give the downward velocity profile 
and the temperature profile at an altitude h ~ h e — s when u u 9i and hi = h a — r are 
known or assumed. 


2.2 The Outflow 

For the radial outflow along the ground the equations for conservation of mass and 
momentum are written 
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where w is the velocity normal to the ground in the z direction and v i3 again the radial 
velocity in the y-direction; e w is an eddy diffusivity for this flow (thus z and z are related 
by z — h 0 — r, and w = — u). 

For the radial flow along ibe ground the effect of buoyancy If neglected and the thermal 
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energy equation is not U3ed. 

The integral form of equation (20) is found by integrating in the z direction. 



or, through the use of equation (10) 
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vydz — constant 

At large distances from the origin a self - similar solution applies cf the form 


v = vi(y/y) p HO 


where £ = z/y 

Substitution into equation (21) gives 


( 22 ) 


Thus p — — 1 and 

„=„,(£) m 

Since, in this analysis, a smooth ground plane is assumed the flow corresponds to that 
of radial fieejet which conserves radial momentum. Again it is assumed that the velocity 
profile 1*(£) has the form 


HO - exp (-c£ J ) 

where c is experimentally determined and is of the order 100. The half width of this radial 
jet i3 


£ - 
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and the velocity profile at large radial distances 13 

u=yi 09 cx? ( 23 ) 

/ 

where yi and v t are reference values to be determined. 

2.3 The Stagnation Region 

In the stagnation region the velocity in the downward flow impinges on the ground 
and is turned radially outward. This outward flow reaches a peak velocity and then 
subsequently decays due to turbulent diffusion. The properties of the downward and of the 
radial outflow are related through conservation of mass and momentum in the stagnation 
region. For simplicity, it is assumed that the radial outflow along the ground (s = o, or 
s — h 0 ) has the form: 


v = 



•for y <y x 


v = Vt 



for y > y 


and that the pressure difference p~ Pea has the form 


(24) 


P~Pc o = ~ v 2 ) for y < y! 

p - Poo = o for y > yi 

Equating the momentum change of the downflow with the overpressure on the ground 
then gives a relation between t>i and y : as follows 
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and substitution for u from equation (17a) and v from (24) gives 
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Similarly equating the mass of the downflow at x = k a with that of the radial flow at 


y = y, gives 




•and substitution for u from equation (17a), and y from (23) gives 


u,, ‘ ; (‘ ■ sj) r c '“' !d < - 'hvi^c^'dt 




Equations (26) and (28) than give 


h 0 \2;ra/ ’ ui \ c ) A.y 

Finally, it is convenient to express tht*3e quantities in terms of references quant Hi-, 
at the height Aj, and to eliminate the virtual origin k 0 . This may be done by using the 
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The expressions in equation (29) can now be written. 
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Also, from (29) 
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8 log 2 A x «i 

la the foregoing expressions b and c are known constants and 9i,D lt «x and h t are 
parameters which describe the microburst at the height A*. 
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(c) At Altitude 

At an altitude hi at a large distance from the grcund, the velocity along the axis of 
the microburst is «x and the velocity profile, ie its variation with y, is given by equation 
(17a): 


-36 / y \* 
u — «x exp — — \^~j 

= Ux exp 4 log 2 Qj-) | 


At any altitude A, the velocity is given by 
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For example, with b — 100, 0 t = I0~ i } Di = 10 s ft, g — 32 ft/sec 2 , u* = 50 ft/sec, 
hi = 10,000 ft then Q — 7. For h = 5,000 ft substitution into (35) gives 


u 

— = .84 exp 



The velocity at the center of the microburst decays rather slowly (ie a reduction of about 

16% at k = 5, 000 ft. from its value at 10,000 ft) and a spreading of the profile radially. 

An aircraft fly : ng through th:3 velocity profile would experience a downdraft and a 

consequent reduction of lift since the effective angle of attach would be reduced by Ac* 'y 

v — . For a microburst of small dimensions (ie approaching the span of the aircraft) 

aircraft 

the aircraft would experience a rolling moment proportional to 77-7^ x where Di 

aircraft 1 

is the diameter of the downburst, in addition to a reduction in lift. Rolling would be in 
the direction toward the center of the downburst. 

(b) Near the ground 

At a distance h above the ground the radial velocity is given by 
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= “(2 log 2) 
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for 



y < /' c V 1+ g 

tfi \8rlos2j Q 
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For example for an aircraft on a path passing through the center of the microburat 
using the same numerical values as in the previous example but with h ~ i.OCC ft., and c 
= SO equations (36a and b) give 
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This velocity profile, shows a rapid decay with altitude above the ground (compared 
with the limiting case k = o). Thus, in thio example, the largest effect of the radial 
outflow is felt at altitudes very close to the ground (ie much less than 1,C00 ft). An 
aircraft descending toward a landing would experience an increase in headwind if moving 
toward the center of the microburst followed by a decreasing headwind close to the center 
and, after crossing the center, an increasing tailwind and an eventual decrease in tailwind. 
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4. Concluding Rcm&xka 

The microburet model described In this paper provides very simple closed form ex- 
pressions for the velocity in the microburst at altitude and in the vicinity of the ground. 
Two constants of the velocity profiles associated with turbulent diffusion may be chosen 
to nt experimental data if these arc available, or used to change the magnitude of velocity 
gradients in the microburet at altitude and near the ground. 

The model, however, does not account for the presence of thermal sources which may 
add energy to the microburst and therefore it may underestimate the magnitude of the 
velocity within the microburst as it approaches the ground. For the purpose of gener- 
ating velocity profiles for flight simulation this deficiency in the model can be offset by 
increasing the effective initial velocity on the axis of the microburet. However, an improved 
model, taking into account thermal sources, is desirable in order to provide a more realistic 
representation of the physical phenomena. 
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